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ABSTRACT 


A flow-type isothermal calorimeter was designed, 
constructed, and tested to directly measure the excess 
enthalpy of binary gaseous mixtures. 

The calorimeter was tested with hydrogen-nitrogen 
mixtures at temperatures of 20, 7, -3, and -73°C and for 
the pressure range of up to 130 atm. The experimental 
data are in good agreement with the literature to within 
the experimental uncertainty, 5 percent. 

The experiment was performed for the nitrogen-carbon 
dioxide system at 40°c and for the pressures up to 130 atm. 
The excess enthalpy data of this system are compared with 
values predicted using the B-W-R equation of state, the 
virial equation of state, the generalized method, and the 
liquid theories. 

The behavior of the excess enthalpy as a function of 
density obeys binary collision theory in the region of low 
densities. In the neighbourhood of the critical point of 
one component of the mixture, the excess enthalpy follows 


the behavior of the enthalpy of this component. 
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NOMENCLATURE 


Constants in B-W-R equation of state 

Constants in B-W-R equation of state 
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Second virial coefficient in volume expansion 
Second virial coefficient in pressure expansion 
Third virial coefficient in volume expansion 
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Pressure drop 
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Gas constant 
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H Constant enthalpy 

sl Component in mixture 

m Mixture property 
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R Reference gas property 
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INTRODUCTION 


Information regarding the properties of materials is 
essential to planning and building which are the busi- 
ness of engineers. A large number of such data have been 
reported from the experimental and/or theoretical studies. 
But the rapid advance of technology into new fields seems 
always to maintain a gap between the demand for basic data 
and its experimental determination. 

For the rational design of processes involving the 
transfer of heat and work accurate enthalpy data are neces- 
sary. Accurate values of enthalpy are directly useful in 
design but they are also needed to test theories of fluids, 
to improve methods of predictions, and to derive other 
thermodynamic properties. For certain fluids and their 
mixtures some enthalpy data which were obtained principally 
from P-V-T data, either directly using thermodynamic rela- 
tions or indirectly using an equation of state, are available 
but those data are limited to a certain region. Particularly 
for the excess thermodynamic properties very few data are 
available. 

For this reason, it is desirable to have excess enthal- 
pies which were measured directly. It is, of course, impos- 
Sible to obtain data for all the materials of interest over 
wide ranges of pressure, temperature, and composition; how- 


ever, experimental data of selected materials and mixtures 


are useful not only in themselves, but also for derivation 
of other thermodynamic properties. This work was aimed to 
obtain data on the isothermal effect of pressure and com- 
position on the excess enthalpy for the mixtures of hydrogen- 
nitrogen and nitrogen-carbon dioxide, and to analyze those 
results for the measurements of the effect of temperature, 
pressure, and composition upon excess enthalpy. 

At the onset of this work, very limited and meager 
information was available on both experimental and 
theoretical aspects of the excess enthalpy of gaseous sys- 


ed: 


tems. Recently information has becomes available on 


both of these aspects through a few ‘lahormategreieg 7 oe eet 


SO. 

The scientific importance of the excess enthalpy is 
that it presents a stringent test of the predictive ability 
of an equation of state and of theories based on intermole- 
cular forces. 

The excess enthalpy, or heat of mixing, is defined 


at constant temperature and pressure, as: 


H = H - 5 x. H.. (a8) 


Figure 1 graphically shows that all the excess enthalpy of 
a binary system is the vertical distance between the 
enthalpy of a real mixture and that of an ideal mixture of 


the same composition. 


ENTHALPY 


Figure. 


~T=Constant 
P = Constant 
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Definition of the Excess Enthalpy for 
a Binary System 
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Because of high accuracies required in obtaining 
the total properties of the mixture and its components, 
direct measurement of property variation is superior to 
investigation of total properties for studying mixing 
phenomena. The total quantities, such as the molar 
enthalpy of a mixture, are also measured by experimental 
procedures, from the P-V-T data, for example. But these 
total values differ, in most cases, only slightly from the 
value obtained by averaging the properties of the mixture 
components. The measurements of the excess enthalpy will 
yield directly the enthalpy deviation of a multicomponent 
system from the ideal mixture case. 

In chapter I, the basic thermodynamic relations 
concerned with this work, a review of the recent litera- 
ture on excess enthalpy of gaseous systems, and methods of 
prediction of enthalpy are presented. 

Chapter II describes the experimental equipment 
construction, measuring instruments used, and method of 
excess enthalpy measurements. 

Comparisons of the experimental results with the 
literature for hydrogen-nitrogen system are considered in 
section A of chapter III. The virial equation of state, 
and corresponding states, such as correlations of Yen and 
Alexander °° ano guid theories>* are tested for their fit- 
ness for nitrogen-carbon dioxide systems in section B of 


the same chapter. Some analytic considerations are also 
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given in both sections of chapter III along with the 


results of this work and those of the literature. 
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CHAPTER I 


PRELIMINARY CONSIDERATIONS 


In this chapter, the thermodynamic relations which 
will be required in this work are presented, a review 
of recent experimental data is given and the methods of 
prediction of the enthalpy of mixtures are considered 
briefly. A list of all experimental work on the excess 


enthalpy of gaseous systems is given. 


ae Thermodynamic Relations 


The first law of thermodynamics, applied to a flow 
calorimeter with negligible potential and kinetic energy 


errecrnis: 


H - H = 2 
Ha ,2, 7 Tala (2) 
x 


where Q is the rate of heat leak, W is the rate of electri- 
cal energy transfer, and F is the mass flow rate. Flow 
calorimeters may be designed for various models of operation. 
In the isobaric mode, the pressure difference, Po-Pi, is 
made small and the fluid is heated to change its tempera- 


ture. Equation (2) becomes: 


F Pp 
W | 29H 
H tS = a oe (x=) dp 
Hn, r,| F oP rT 
pay x 


1! a hy 


(3) 
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where the integral term is the correction for the small 
pressure drop and heat leak is assumed to be negligible. 

In the isothermal mode, a pressure drop is imposed on the 
fluids and electrical energy is added to return the out- 
let temperature to that of inlet. It is possible to 
utilize this scheme only when the Joule-Thomson coeffici- 
ent is positive, i.e., when the fluid cools upon expansion. 


In this case Equation (2) reduced to: 


where Q is assumed negligible and the integral term corrects 

for any mismatch between the inlet and outlet temperatures. 
In a flow calorimeter operated in the isenthalpic 

mode, no energy is added to the system and the heat leak 


is made negligible. For this case Equation (2) reduces to: 


= 


H - H =O) (5) 
Ho, et 


The relations presented above involve integral 
changes in enthalpy. This integral data in a single phase 
region may be differentiated to yield the derivative 
enthalpy properties: the heat capacity, Cp, the Joule- 
Thomson coefficient, u, and the isothermal throttling 


COerEiciert, ¢. Thus: 
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These three derivatives are related by the mathematical 


identity: 
Cae UCD. (9) 


The isothermal throttling coefficient, ¢, can be 
expressed in terms of volume, temperature, and pressure 


as: 
atyoe pe (ey. 
a on (10) 


thus allowing this coefficient to be obtained from compres- 
sibility data or an equation of state. 
The virial expansion which is an open form of an 


equation of state, is a power series in density: 


a ae (10) 


“ 
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‘ 


Thewrerms Bb, Coe... are known, as’ the second, third, ... 
virial coefficients and they are related by the statis- 
tical theory of imperfect gases to the interactions of 
molecules in pairs, triples,.etc. 

A similar power series is an expansion in pressure: 


PV 
RT 


1 +) BPs |] ic pe Ci) 


The coefficients of two series are related by: 


B! = B/RT (133) 
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(C-B PRE (14) 


For sufficiently low densities, the virial equation of 


state truncated after the second term can be expressed 


as: 
PV _ 
ar OUT l e+ B/Ve (lla) 
Equation (lla) can be applied to a binary mixture and 
its components. From the definition of thermodynamic excess 
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and with the thermodynamic relations: 


9H aV 
(—=) = V - T (==) (10) 
oP oT, 
and 
He SS Pieri = H_) (1a) 
m INIA 5B a 


one can show that 


Ho ah - | dP Gin) 


ae oT 
and 
P E 
de = [fat om Ged an (18) 
oT 
‘e) Pp 
E ap™ 
Moreover, since B- and ap are a function of tempera- 


ture only, we have: 
E 
Eau Loe dB 
H —e" E Al (om | P. (19) 


The mixture virial coefficient, oa is composed of the 


contributions of three kinds of molecular interactions: 


a Zz 2 
Bn = Xn Ban + 2xX,XpBap + Xp Bap GeO) 


where the coefficient BAR 1S,.due, to. the interaction. of 


molecules A and B, and Ban and Bap are due to A-A and 


B-B interactions, respectively, 
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Substituting Equation (20) into (16), we have: 


pairs 
B = 2X,XpE CU2:1%) 


where E, the excess second virial coefficient, is de- 


fined as: 


E= Be - F i mpited yor (22) 


Thus for low densities, since E is a function of tempera- 
ture only, Equation (21) can. be substituted into Equation 


(19) giving, for the excess enthalpy: 


Fi) a genes PIE = 2 (Gp ile (23) 


The enthalpy change on mixing can be determined 
in a flow caiorimeter by mixing two pure gases in a 
mixing chamber and adding electrical energy to equalize 
to inlet and outlet temperatures. The first law of 


thermodynamics for this case may be expressed as: 


m = oe leh Q 
: i. a = ° (24) 
Tap, By ei THe ADH 


Correction may be made for the differences in pressure 


1 
= 


Hy 


and temperature between the inlet and outlet. For the 
isobaric and the isothermal mode of Equation (24), the 


excess enthalpy can now be determined by: 


~ 


a0 
E E “| 
H = H ii aan | (5) 
am i "5 prix’ 


As shown in Equation (1), needed in the calculation 
of the excess enthalpy are enthalpies of the pure components 


and that “oftheir mrxture. 


B. Experimental Enthalpy Determination 


The enthalpy has been obtained experimentally through 
the measurements of Cp, u, ¢ as described above. Masia? 
presented a survey of these various methods for collecting 
enthalpy data, and a review on the Joule-Thomson effect was 
covered by Johnston?° , Potter “® has surveyed both Joule- 


Thomson and isothermal throttling coefficient data. 


C. Enthalpy Determination Through Empirical Methods 


Since it is inconceivable that information on the 
physical and thermodynamic properties of all known substan- 
ces and their mixtures will ever become available, it is 
essential to have methods to estimate and predict values for 
desired properties. In this regard, it is necessary that 
the method be subjected to and defended against reliable 
experimental results. 

The literature presents numerous methods for esti- 
mating enthalpies of pure substances and their mixtures. A 
comprehensive review of prediction methods is given by 


ZO 2 


Hobson and Weber, and recommendations on the various 
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correlations are discussed by Reid and Sherwood>*. 


Lydersen, Greenkorn, and Hougen>° have presented 
generalized charts based on critical properties of gases 
for estimating thermodynamic properties. 


Various equations of state, such as the Benedict- 


Webb-Rubin (eoweoR ela. Redlich-Kwong meee. Hirschfelder- 
Buehler-McGee-Sutton TS ES APRS and the virial equations, 
can be used to determine enthalpies. A recent review for 


B-W-R, R-K, and H-B-M-S equations of state has been present- 
ed by Yorizane and Masuoka’?. Recently Senter. has developed 
a method based on the principle of corresponding states using 


molecular parameters to predict properties of a system. 


D. Prediction Methods of Mixture Enthalpy 


For gas mixtures, deviations in behaviour from an 
ideal gas mixture may be caused by either or both of two 
effects: 

1) the component behaviour may be nonideal, and 

2) the interactions between components may contri- 

bute substantially to the nonideality. 

It is possible for a mixture to be nonideal even 
though the behaviour of the individual pure components 
is ideal at the same temperature and pressure. 

In most engineering problems, enthalpies must be 
generated from correlations and prediction methods. The 


experimental thermodynamic properties of gas mixtures have, 
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until recently, followed the same line as that taken for 
pure gases. The number of systems of interest is so 

large and the rate of production of experimental enthalpy 
data so small because of the time and effort involved 

that it is unlikely that experimental data are available 
for the problem under consideration. The ideal solution 
would be the calculation of thermodynamic properties, in- 
cluding enthalpy, from a knowledge of the interactions 
between molecules. At present this can only be done for 
relatively simple molecular models, but increasing efforts 
are being made in this area. For instance, Hermsen and 
Peeusnite-— and Eckert, Renon, and Pransnier. have calcu- 
lated excess functions for binary liquid mixtures of 
hydrocarbons at low pressure which are in agreement with 
experimental values. Osborne ® has also presented a method 
for the prediction of liquid mixture enthalpies from a mole- 
Culan modelsitor liquid mixing. 

When attempting to determine the enthalpy of a mix- 
ture by means of a semiempirical method, one seldom has 
available the mixture parameters required to directly apply 
the model. For this reason, a prescription is required 
to obtain the parameters of the mixture from those of the 


pure materials. 
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Recently Nathan@4 has reviewed the types of mixture pro- 


perty prediction methods. 


1. Pure Components 

The assumption of zero heat of mixing allows the 
enthalpy of a mixture to be calculated from the sum of 
the enthalpies of the pure components present. This 
method may give good results at moderate pressures and 
relatively high temperatures, for chemically similar com- 


pounds such as hydrocarbons in a homologous series. 


2. Equivalent Component Concept for Mixtures 


Several methods have been proposed which consider 
the mixture to be a hypothetical pure component and corre- 
late enthalpies based on a parameter such as model aver- 
age boiling point, mixture molecular weight, or mixture 
specific gravity. 

The assumption is made that a mixture and a pure 
component with the same correlative parameter exhibit the 
same enthalpy behaviour. Canjar and Peterkal® have pre- 


sented correlations of this type for hydrocarbons. 


Buy Corresponding States 


The principle of corresponding states was first 
applied to the correlation of P-V-T data, but the extension 
to enthalpies using the integrated form of Equation (10) 
follows directly. Among the earliest correlation of this 


types are those of Cope, Lewis, and Weber? + and Edmister??, 
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Later work on generalized correlations has employed a 
third parameter to improve the agreement with experi- 
mental data. Lydersen, Greenkorn, and Hougen~° used the 
compressibility factor at the critical point, Zc, while 
Gus land Pitzer!> employed the acentric factor, w, which 
is related to the shape of the reduced vapor pressure 
curve. In applying these correlations to mixtures it is 
necessary to obtain values for the critical temperature, 
Tc, the critical pressure, Pag, and the third parameter of 
the mixture. In most cases, recourse is made to rules 
which allow the mixture critical properties to be calcula- 
ted from those of pure components. For example, the 


values for a number of correlations are obtained by the 
Sisk 


linear mixing rules suggested by Kay, but nonlinear 
mixing rules have been proposed by Pitzer and Hultgren*’, 
and Prausnitz and Gunn??, 

4, Enthalpies from Fugacities 


The effect of temperature on the fugacity of a com- 
ponent in a mixture is given by: 


—————— 2) 
aT P,x RT 


and if relation for the fugacities is known, the partial 
molal enthalpy H, can be calculated. The disadvantage of 
this method is that any error in the expression for 


fugacity, which is derived from P-V-T and vapor-liquid 
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equilibrium data, is magnified in the differentiation. 
The advantage claimed for this method is that the enthal- 
pies obtained are consistent with the vapor-liguid equi- 
librium data, a point of importance in multicomponent 
mixtures. Klein? has used this method inversely to cal- 


culate fugacity coefficients of methane and nitrogen from 


their mixture enthalpies. 


or Equation of State 


The enthalpy of many substances is known in the 
ideal gas state from spectroscopic calculations and direct 
experimental determinations at low pressures. With the 
effect of temperature on enthalpy known at zero pressure, 
the effect of pressure on enthalpy can be obtained from 
Equation (10) if an equation of state is used to represent 
the P-V-T behaviour of the substance. 

A large number of equations of state have been pro- 
posed, but only a few have been used for extensive calcu- 
lations of enthalpies. Martin>© has recently reviewed 
equations of state primarily from the standpoint of repre- 
sentation of the P-V-T surface of a substance. Recently 
Single equations like the Redlich-Kwong>~ equation have 
been used for enthalpy prediction and the recent modifica- 


ager yey) 


tion by Wilson is claimed to improve the accuracy of 


such predictions. More complex equations such as the 
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Benedict-Webb-Rubin” and Martin-Hou equations have also 
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been used to predict thermodynamic properties including 
enthalpy. 

A problem arises in the extension of equation of 
state calculation to mixtures since the constants for 
particular mixtures are not available and usually must be 
calculated from empirical mixing rules utilizing the con- 
stants of the pure components. The virial equation is 
one of the few equations for which the composition depend- 
ence of the constants is known exactly. This equation is 
restricted to vapors since the power series relation tends 


to diverge at densities approaching those of liquids. 


E. Importance of the Excess Enthalpy 


A knowledge of the precise thermal data is required 
for the accurate design of numerous types of process equip- 
ment and. the.lack of sufficiently sound data forces the 
design engineer to overdesign with a resultant increase in 
costs. The enthalpies for pure components are usually 
known to a good degree of accuracy, but the same is not 
generally werue for mixtures. g Accuracyean the excess enthal- 
py permits economical operation in gas separation, and 
gas purification, and cryogenic process by improved design 
methods. 

Accurate data are difficult to determine for a multicom- 


ponent system. Even for the simplest and most common systems, 
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multicomponent data for fluids are almost nonexistent. 
Engineers and scientists required to employ such data 
usually must resort to simple mixture rules or correla- 
tions developed for specific situations to predict mixture 
methods of approximation. The validity of use of those 
methods is also in question when applied to the systems 
other than the simplest nonpolar molecules in the regions 
of moderate high density. 

It appears then that, in the field of gaseous mix- 
tures, there is a definite need for thermodynamic proper- 
eVec ume Data or hic type mot only aresot direct utility 
to the process industry but also can be further used to 
obtain valuable information of a more scientific nature. 

One of the classic problems of both physics and 
chemistry has been that of understanding the molecular 
forces and their effect on gases and liquids. Because of 
the fact that the excess enthalpy reflects directly all 
of the imperfections of a solution from that of the ideal 
case, it is a very useful quantity for testing models 


based on the intermolecular POrCeson. an 


F. Necessity of Direct Determination of Excess 
Enthalpy 


Deriving the excess enthalpy from enthalpy data 
obtained through experimental methods of pure components 
and mixtures requires taking the differences of the total 


quantities, Ha and DX sH ine Equa tons), which ~ane. ip, 
a 


cm 


general of the same magnitude. 
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When the subtraction is 


carried out to obtain ce the uncertainty in the excess 


enthalpy is the sum of the uncertainties in the individual 


enthalpies. Since the 


Therefore, 
compared to the individual values 
duced in obtaining H® in this way 
percentage of the excess enthalpy 


enthalpies of pure components and 


excess enthalpy is small 
used, the error intro- 
is usually a significant 
LE 


value calculated. 


their enthalpy change on 


mixing are available, the enthalpy of mixture is obtained 


from Equation (1), as: 


Hie Hae 
AL 4. 


G. Background of This Work 


The gaseous systems, which were investigated in this 
work are hydrogen-nitrogen system and that of nitrogen- 


carbon dioxide. These systems are ideally suited for an 


exploratory investigation of this nature because, from an 


economic viewpoint, they are cheap materials, and froma 


molecular point of view, since their molecules are nonpolar 


and simple, it fits the assumptions involved in the Lennard- 
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Jones potential quite well ~. Pure component enthalpy data 


pale Si) 


on both hydrogen” and nitrogen are in abundance and 


enthalpy data are available on at least some compositions 


of their mixtures. A number of experimental data for the 


thermodynamic properties of carbon dioxide are also available 
9,45 


easily - Moreover since excess enthalpy data experimentally 


I= 


Zhai, 


determined for the hydrogen-nitrogen system are available 
from Knoester, et Sey, this system is used to test the 
equipment constructed for this work. 


A complete review of experimental thermal data on 


various pure substance and mixtures was given by Mather +, 
and Yesavage, et alae Knoester, et Piles presented a 


compilation of thermodynamic properties, both experimental 
and calculated, of simple molecular gaseous mixtures, and 
Van Eijnsbergen and Beenakker>° listed tabular and graphical 
data of thermodynamic properties for cryogenic fluids such 
as normal hydrogen, para hydrogen, helium, neon, nitrogen, 
carbon monoxide, oxygen and argon as well as methane and 
xenon. 

Table I shows a review of all the experimental excess 
enthalpy data up to date for the gaseous mixtures. Theoreti- 
cal determination and indirect measurement of the excess 


enthalpy from the other excess thermodynamic property by 


experiment were also published as listed in Table II. 


it 
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Table I 


Experimental Excess Enthalpy Data for Gaseous Systems 


Year System Authors References 

1962 Hydrogen-nitrogen Beenakker and alk 
Coremans 

1965 Hydrogen-nitrogen, Beenakker, et al. 2 


hydrogen-argon, 
nitrogen-argon, 
methane-hydrogen, 
argon-methane 


1967 Hydrogen-nitrogen Knoester, et al. 35 
hydrogen-argon, 
nitrogen-argon, 
hydrogen-nitrogen- 


argon 
1968 Methane-argon, Van Eijnsbergen and ays) 
methane-nitrogen, Beenakker 


hydrogen-methane, 
helium-methane, 
helium-argon 


1969 Nitrogen-methane Klein 32 
1969 Nitrogen-methane Klein, et al. 38 
1970. Nitrogen—carbon Hejmadi, et al. 20 
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Calculated Excess Enthalpy Data for Gaseous Systems 


Year 


TENS Gi) 


1968 


19:69 


Table II 


system 


Hydrogen-nitrogen, 
hydrogen-argon, 
nitrogen-argon 


Methane-argon, 
methane-nitrogen, 
hydrogen-methane, 
helium-methane, 
helium-argon 


Hydrogen-nitrogen 


Authors 


Knoester, et al. 


Van Eijnsbergen 
and Beenakker 


Hsi and Lu 


References 
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CHAPTER II 


EXPERIMENTAL METHOD 


In this chapter, the equipment construction, 
measuring and controlling instruments, and the method 
of measurement for the direct determination of excess 


enthalpy are described. 


A. General Survey 


The principle upon which the direct determination 
of the excess enthalpy based is the fact that when two 
different gases which are initially at the same temperature 
and pressure are mixed, a temperature change will occur. 
Since most gaseous systems have a positive excess enthalpy, 
heat can be supplied to the system to compensate for the 
cooling effect due to mixing. Hence an isothermal flow- 
type calorimeter is suited for these studies. 

Figure 2 represents a schematic experimental method 
for a continuous flow system. Two pure gases A and B 
which are maintained at the desired pressure, pass through 
the heat exchanger where their temperatures are equalized 
and adjusted to that desired for the run. They then enter 
the mixing chamber. The inlet temperature, before mixing, 
and that of outlet, after mixing, are measured and balanced 
by venergy input. | The flow rate of the mixture, the compo- 


sition, and the rate of input energy are determined. The 
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heat of mixing is then calculated from the quotient of 


the input energy and the flow rate. 


B. Experimental Equipment Construction 


The apparatus is consisted of four main parts, 

namely, 

(1) a calorimeter, 

(2) a low temperature bath, 

(3) a room temperature bath, and 

(4) a gasometer. 

The detailed schematic flow diagram of the appara- 

tus is shown in Figure 3. Details of the construction 
are as follows. 

l-)) Galorimeter 

As shown in Figure 4, the calorimeter which is set 

in a stainless steel vacuum jacket, has a mixing chamber 
(The detailed drawing for this is shown in Figure 27 in 
Appendix A) consisting of: 

(a) three concentric inner copper shells to accom- 
plish complete mixing, 

(b) one outer stainless steel shield to minimize 
heat exchange between the mixing chamber and 
the surroundings by radiation, and 

(c) a 16-foot-long 23 gauge nichrome wire which is 
approximately 20 ohms at room temperature, was 
wrapped over the first and second copper shells 


for heat supply to the gas mixture. 
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LEGEND 


1, A gas Container(s) z 
2. B gas Container(s) 

3, Shut-off Valves for Pure Gases 
Constant Low Temperature Bath 

5. Thermocouple Vacuum Gauges 

Vacuum Valves 

Vacuum Release Valve 

Air Cold frap 

Vacuum Punp 

Calorimeter 

Controlled Heat Input 

Liquid Nitrogen Cooler 

Liquid Nitrogen Dewar (50 1,) 
Pressure Reducer 

High Pressure Mercury Manometer 
Stirrers 

Sensors for Heat Input Controller 
Platinum Resistance Thermometer 


Figure 3. Schematic 


Room Temperature Water Bath 

Nitrogen Vapor Exhaust Control Valve 
Shut-off Valve for Mixture Stream 
Mixture Flow Rate Regulating Valve 
Sampling Valve 

Flowmeter 

S-way Solenoid Valve ; - 

Release Valve for Mixture Exhaust 
Water Manometer 7 

Thermocouple for Measuring Mixture Temperature 
Check Valve 

Balancing Weight 

Water Drain Valve 

Water Level Indicator in the Gasometer 
Gasometer 

Heise Bourdon Tube Gauge 

Switching Valves for Pressure Reading 


Flow Diagram of the Apparatus 
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Schematic Diagram of the Calorimeter 


Figure 4 
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Three 6-junction copper-constantan thermopiles were con- 
structed and mounted on the calorimeter as shown in 
Figure 4. For this thermopile construction 40 gauge 
copper wire and 32 gauge constantan wire were used, and 
the distance between hot and cold side terminals was 
made 1% feet which is long enough for negligible small 
heat conduction from the hot side to the cold side through 
thermocouple wires. The thermopile sheath was made of 
4-inch-long and 3/16-inch outside diameter thin wall stain- 
less steel tubing. The function of each pair of these 
thermopiles in Figure 4 is as follows: 
pair A-D: to check if inlet gases to the calorimeter 
have reached the same temperature as the 
controlled bath temperature, 
pair A-B: to monitor the energy supply, and 
pair B-C: to check if mixing is completed in the 
mixing chamber. 
For the quick response to the temperature changes 
of the gas stream to the thermopile junctions Apiezon T 
was filled into the gap between thermopile sheath and 
the junctions. A 3%-foot-long 3/8-inch outside diameter 
stainless steel tubing was spiraled and connected by weld- 
ing to the outlet of the mixing chamber and the thermopile 
pair B-C was mounted at both ends of this spiraled tubing. 
A Conax gland was used to bring out the copper lead 


wires of the heater from pressurized mixing chamber to the 


es 
-nOo Sisw eeliqomielt asinegenoo-19eqqo> .noitonut+o sesat 
ak nworle es sustomizols> sig no bodavom bas besoumte 
eids 10oF .& oteplt 


spyne Ob noLtoust: 
- Dns stiw teqgoo 


,DSsau siIsw stiw netces2aimoo SPUBP > 
od meawted somsatetb ed? 


oo, siigqonirsns 


ee 


Bits 
esw. alaninoyves obi? FLoS Ons 2 
{isme efdipil pan 193. dpuons pnof ef doidw tee% ¢l ebsm 
dpwowts shie blood sdx ot sbhie tod. ods.noxd noissubaes seed 
70 Sham aw dtserie sltqomusd? sit 2oxtw silquovemzeds 
-ntste Lisw cairo fogi-dl\t Bas porol<domi—-b 
sasgns ik tsq dose to nott+snvt oAT .pnidet Ieaete aged 
wolflot 2 vad bh suupldi at esligommeds 
tevemiiolss sit oF esesp asink t+ Aasdo of i~A TIO 
see Sit bserloss1 sven 


Srussrsame3 


w361i9qmSes At6sd Nellortnoo 


Sit tos f£n0m od 
assdi>d od +O-H tis 


i+ pe 


‘ 
.viqque yptens :G@-A xt8g 


wits 
MHLixXoiM If 


sit mi Hesalqmon ear 
paLxrm 


asdmsoo 
Qi senoqeoy Aoiup eit 169 
ot mesit2 esp eft Yo 


Bboliit esw 


9 ry ad 
a3 


BOpnedio siutsreqme: 
? noseigA encisonavt sligomisds 9 
SifgGomisnis nsswied gsp sd3 osnt 
domi-—8\E pnol-toot-#1£ A -enotsonavt sds 
& Belsitge asw pnidus foesca gasinisses 


to toftuo ert of pat 
Setavom asw O8@ akeg 


bas Assesde 
UO 


ta2enisrb sbiss 


7 


“Disw yd bescsrinos 
pNIxim edt 


tsdmeno 


9f sgomtedt sit Bris 
spnides Oslarige eida ito abnas ritod 45 
Bbasip xsAod A 


beel teqgom oft tuo poirzd of Bsey asw 
edt oF xedmedo paixim besixueeexq mort tsdssd sid Fo aeriw 


30 


evacuated vacuum jacket. To bring the lead wires of the 
thermocouples and the heater from evacuated vacuum 
jacket to the atmosphere, an 8-pole vacuum seal was 
mounted on the 2-foot-long %-inch outside diameter stain- 
less steel tubing which was welded on the lid of the 
calorimeter vacuum jacket. This jacketwas installed for 
minimizing convective and radiative heat exchange from or 
to the mixing chamber with surroundings. The jacket was 
constructed with a piece of stainless steel pipe of 9-inch- 
high and 4%-inch outside diameter. For fixing the lid, 
a Teflon O-ring was used with 16 stainless steel bolts to 
tighten. 

2. Low Temperature Bath 

It is necessary to maintain a constant temperature 
level during the experiment. For this purpose a stainless 
steel Dewar was constructed and insulated with styroform. 
These tmcapactly, Ofecthe baths aboutiG/ liters. To intro—- 
duce both inlet gases to the calorimeter with the equal 
temperature, 50-foot-long soft copper coils were installed. 
Beside the calorimeter Jacket a platinum resistance ther- 
mometer was located to measure the absolute temperature 
of the bath. The temperature controller used a thermistor 
probe as a sensor and heaters of either 300 or 500 watts. 
The bath was well-stirred by a 2-stage impeller. 

3. Room Temperature Bath 

The temperature of gas mixture was required to be 


at the room temperature before it reached the regulating 
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valve and the gasometer. A bath of about 80 liters was 


installed for this end. 


Cc, Measuring and Controlling Instruments 


The major instruments to measure and control various 
values to determine the excess enthalpies are described 
briefly. 

(1) The absolute temperature in the low temperature bath 
was measured by a platinum resistance thermometer and 
a K-5 potentiometer which were calibrated by the U.S. 
National Bureau of Standards. Their calibration data 
are presented in Tables VIII and IX in Appendix B. 

(2) The pressure of the system was measured by Heise 
bourdon tube gauges. No difference between the dead 
weight gauge and the Heise gauges were found through- 
out the entire pressure range. 

(3) A thermocouple vacuum gauge was used to read the 
vacuums in the stainless steel Dewar and the calori- 
meter jacket. The vacuums were as good as 2 microns. 

(4) The electric energy was supplied to the gas mixture 
by a low voltage d-c power supply. 

(5) The input energy was measured by a 2-channel strip 
chart recorder which is in an error of less than 0.5 
percent in full scale (including reading error). 

(6) The required time to collect the mixture in the gaso- 
meter wasS measured by an electric timer with an accu- 


Lacysot 0.01 minute. 
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(7) To read the null point of the thermopile pairs a 
d-c null voltmeter was used. 

(8) The pressure drop across the calorimeter was checked 
by a high pressure mercury manometer. 

(9) The bath temperatures were controlled within +0.05°C 
by temperature controllers. 

(10) An 0.1 standard resistor and a 10 ohm standard resis- 
tor were used to construct the platinum resistance 
thermometer circuit and calorimeter heater circuit, 
respectively. 

(11) A 150-liter-gasometer was used to measure the mass 
flow rate of the gas mixture. The calibration data 
for the gasometer are given in Table XIII, Appendix B. 

(12) Six-junction differential thermopiles were used to 
investigate temperature differences. All these 
thermopiles indicated an error of less than +2 
microvolts (equivalent to +0.0083°C) at any tempera- 
ture level of this work. Therefore these were used 


WithouUsesCcalLabrations. 


D. Method of the Excess Enthalpy Measurement 


The experimental determination of the excess enthal- 


py was performed as follows. 
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Firstly the pressures in the both pure gas con- 
tainers were measured. 

Opening the valve of the higher pressure side, the 
system was filled out with that gas. Because of the gas 
compression effect in the system, a temperature gradient 
appeared in the calorimeter. It took about 4 hours to 
reach thermal equilibrium in the calorimeter, but opening 
the regulating valve to allow gas to pass through the 
calorimeter, reduced the time to about one hour. When 
thermal equilibrium is obtained, the reading of d-c null 
voltmeter indicates zero and at the same time this is 
indicated on the strip chart recorder. 

When the pressure of the higher side reduced to 
the same pressure as that of the lower side, the valve of 
the lower pressure side is opened. The two pure gases 
flow into calorimeter and mixed there at same pressure 
erg ucated by the pressure gauge. At the same time a 
temperature difference occurs between inlet and outlet 
of the mixing chamber. This deviation is recorded on the 
recorder after amplification by a d-c null voltmeter. 
Then the required energy is supplied by a manually con- 
trolled d-c power supply, and the voltage drops across 
the heater in the calorimeter and the standard resistor were 
recorded on the two-channel recorder. From time to time 
the bath temperature was checked, and also the temperature 


of the inlet gas to the calorimeter was compared with that 
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of bath fluid. The mixture flow rate was controlled by 
a micro-regulating valve so that there was negligible 
pressure drop across the calorimeter. 

In order to measure the flow rate of the gas 
mixture, the mixture is collected in the gasometer at 
the standard condition, and the time is measured by an 
electric timer. The flowmeter indicated roughly the 
velocity of the mixture stream. A 3-way solenoid valve 
was installed between the flowmeter and the gasometer to 
choose the flowing route of gas stream to the gasometer 
or to the atmosphere for vent. 

All bath temperatures of each bath were controlled 
at their desired temperature levels by each controller. 
Bath fluids were water for above 0°C and 2,2,4-trimethyl- 
pentane for lower temperatures. 

1. Composition Preparation 

The composition of the mixture is determined by the 
ratio of the volumes of both pure containers combined, since 
those are always at the same pressure when emptying slowly 
during the measurements. But even if a set of known net 
volumes of gas cylinders is combined for the run, it is 
difficult to obtain a supply of constant composition mix- 
ture because the compressibilities are not the same for 


the different gases. 
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From the equation of state we have: 


PPV 
n= SRT: (27) 
When gas is released from a constant volume container 
at constant temperature, the change of the number of 


moles in the container is: 


P P 
2 tI V 
Ang=* (—=/- —=y"= (28) 
Zo Za RT 


where Z4 and Zo denote compressibility factors at pres- 
sures Po and Po at same temperature T. 


But for a very small change of pressure, we may 


assume Z) = Zo, then Equation (28) can be expressed as: 
V 
An = AP SRT (29) 


where z is the compressibility factor at T and average 


value of Pi and Po. 


For a binary system formed with gases A and B 
which are supplied to be a mixture as the same pressure 
drop AP, the mole percent A in the mixture can be cal- 
culated from the following equation: 


Vee "2, (P,T,) LOG 
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mole % A(P, Tye 


(30) 
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where Va? net volume of A gas container, 
Vp? net volume of B gas container, 
Ty? absolute temperature of gas A to be 
supplied, 
T,? absolute temperature of gas B to be 
supplied, 


2 Pl) COMP Geese bh] bityeLactorsot. gaseA at 
A A 
Peace cnc, 


Z,(P,T,): compressibility factorfof gas B at 
beend TR: 
2. Temperature Control and Its Measurement 
To have the temperature ievels lower than room 
temperature, liquid nitrogen was supplied into the bottom 
of the %-inch outside diameter copper coil in the bath, 
boiled off after absorbing heat from the surroundings, 
and the gaseous nitrogen was vented to the atmosphere. 
To minimize consumption of coolant and the temperature 
fluctuation, a pressure reducer and a regulating valve 
were used for control of liquid nitrogen evaporating rate 
Prec iemscOUDemECOlL lL Ofethenbath,. 
The absolute temperature of the bath was measured 
by the platinum resistance thermometer with the aid of 
a potentiometer and the conversion table based on the 
Callender equation. The circuit diagram for this measure- 
ment is shown in Figure 5. For use of the conversion 


table the following equation was used. 
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Platinum resistance 
thermometer in low 
temperature bath 


0-1,000 ohm 


adjustable resistor 
(Helipot type 

10 Turns 

+ 0.1% linearity) 


lOohm 
standard 
resistor 


On- off 
switch 


NOTE: 


E,; and E> are connected to the potentiometer. 


Figure 5. Platinum Resistance Thermometer Circuit 
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38 
ER, 
SHE Toe = Y (315) 


where = oh voltage drop in volts across the platinum 
resistance thermometer, 


E.: voltage drop in volts across the 10-ohm 
standard resistor, 


R,: resistance of the 10-ohm standard resistor, 
R_: absolute resistance of platinum resistance 
P thermometer at 0°C (= 25.564 ohms). 


The Y in Equation (31) was compared with the values 
of the conversion table and converted to the temperature 
in degrees centigrade. 

3. Pressure Measurement 

The pressure was read by a Heise bourdon tube gauge 
when thermocouple pair A-B read zero, i.e., the deviation 
record was passed over the zero line. To obtain the exact 
pressure reading, actual performance was as follows: 

When mixing was taking place in the mixing chamber, 
the temperature dropped. Then such an amount of heat was 
supplied such that this cooling was nearly compensated 
(line c in Figure 6). If the heating power is not changed, 
the temperature of the outlet of the mixing chamber will 
slowly rise (line d) because the pressure in the system 
decreases gradually, which means a lower gas flux and al- 
most always a lower molar excess enthalpy. At the inter- 
section of line dy with the zero line the cooling effect 


is compensated (measuring point 1). At this moment the 
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Figure 6. 


Differential Voltage Recording 
for Energy Input to the Gas Mixture 
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set value of power was taken and at the same time the 
pressure was read. A second measurement was performed 

at somewhat lower pressure by diminishing the heat supply 
by a small amount (line e). In this way a set of data 

at various pressures can be obtained in a rather short 
time period. 


4. Flow Rate 


The mass flow rate of the gas mixture was measured 
by collecting gas mixture at a standard state in the 
gasometer. Following equation was established for the 
calculation of flow rate: 

(,,.003,54x Ome XaVe X90 XnE 
m m 


Fo= —_—————__ __ _-.-________—____—-. (32) 
(243% 15+ Ta eZ Sil 


m 
where F: flow rate of gas mixture in gram-moles per second, 
Vint the volume in cubic centimeters for one centi- 


meter gasometer rise, 


as gasometer rising distance in centimeters during 
Q seconds, 


P_: the partial pressure of the mixture in the gaso- 
meter in millimeters mercury, 


T : temperature of the mixture collected in the gaso- 
meter in centigrade, 


z : compressibility factor of gas mixture at Ta and 
Ee 


m- 
The partial pressure of the mixture, Pal is expressed 


ass 


=— ee £.. * P (33) 
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where Pri total pressure in the gasometer, 


Put partial pressure of water vapor in the 
mixture of gasometer, and 


Py: pressure reduction equivalent due to the 
gas mixture diffusion into the water 
through the small contact area in the 
gasometer. 

It was found that there was not a sensible pres- 

sure change due to the Pu and Ph: Therefore Equation 


(33) becomes: 


P = .P = 


m sit Pharometric pressure CE 


To make use of Equation (32), ve =e2yiore ece7em (refer to 
Table XIII in Appendix B) and d = 63.1 cm were chosen as 
constants. Ta and Pa were room temperature and barometric 


pressure, respectively. 


5. Input Energy to the Gas Mixture 


Heat dissipated across the nichrome wire in the 
calorimeter is measured by the circuit as shown in 
Figures) .. lLhescurrent of the circuit can. be written. ase 

I= °2 pp (35) 

Owl 
where I is current in amperes, E, voltage drop in volts 
across 0.1 ohm standard resistor (Ro 4 = 0.1 ohm). 


Then the power, W, in Joules per second, is expres- 


sed as: 


(36) 
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Figure 7. Calorimeter Heater pac Wal ae 
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where Ey is voltage drop in volts across the heater in 


the calorimeter. 


Combination of Equations (35) and (36) gives: 


Wae=) (10)-x E, x E, - (370 

Since the excess enthalpy, Hee Ln Joules per gram-mole, 
is calculated from the quotient of the input energy and 
the flow rate, or from Equations (24) when neglected 


heat leak, we have: 


= 


Hy] = 


(38) 


. i 
row / 7 eae" - or _ 
a 7 

aa © _~ - w Te 7 j 
4 gt aetsed oft 2a0%06 adtov ni qoxb epeefoy mk F alee 


_sedemixoin> 8 


cmmek (9€) Bas (dé) enoitsupa to noitsnatdmod 


(ve) «x fx (OL) = w 


,alom-msxzp isa eelvol ai a .vgisdine e2eaoxs edt eonte 


= 


bas ypiens sugai sx to gasistoup sit moxr2 bet siuolae , 


bessefpen aedw (NS) enoitsupa moxt xo \sdex wolt exis 
{ ‘ 
‘oven ow ,xeeL a: 


a. a 


(8£) 


44 
CHAPTER iff 


EXPERIMENTAL RESULTS AND THEIR ANALYSIS 


This chapter covers the experimental results of 
this work, comparisons with the literature, and data 
analysis. 

Presented in section A for hydrogen-nitrogen sys- 
tem are the smoothed excess enthalpy data, comparisons 
with the literature, and their interpretation. 

In section B experimental results for nitrogen-carbon 
dioxide system are reported, and comparisons are made with 
the virial equation of state truncated after the second 
term, the B-W-R equation of state and the corresponding 
state correlations of Yen and Alexander®? and liquid 


wHeoriesos4 


A. Hydrogen-Nitrogen System 


The direct determination of the excess enthalpy of 
the hydrogen-nitrogen system has been covered over the 
pressure range of 5 to 130 atm., the temperatures of nomi- 
Nale20;,. /;,2—-3,. and Safes and for the compositions of 
nominal 26, 51, and 76 mole percent nitrogen. The experi- 
mental raw data of the excess enthalpy as a function of 
pressure and composition, are presented in Table XVII in 
Appendix D. A detailed example for the calculation of the 
excess enthalpy from the experimental measurements is given 


in Appendix E. 
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lL. Material Used 
The sources and purity of the gases used are given 
in Table III. The gases were used without further puri- 


frcativon-. 


Pes Experimental Work 


The experimental method was followed as the same 
way described in section D of previous chapter. 

The pressure of the hydrogen containers was between 
1,800 and 2,000 pounds per square inch and that of nitrogen 
containers between 2,100 and 2,300. Both pure gas containers 
were at room temperature. They were regular standard com- 
mercial gas cylinders HK (calculated as 1.51 cubic feet in 
net volume) for hydrogen and K (calculated as 1.52 cubic 


feet in net volume) for nitrogen. 


Composition Preparation 


For the hydrogen-nitrogen system, since for both pure 
gases their compressibility changes with pressure at constant 
temperature are almost linear, the composition was deter- 
mined by calculation method. 

Since both pure gas cylinders are set in the same 


room temperature, Equation (30) reduces to: 
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Table III 


Sources and Purity of the Gases Used for the 
Hydrogen-Nitrogen System 


(1) Nitrogen (High Purity) 
Supplier: Consumer's Welding Co. 


Typical Analysis Data: 


nitrogen 99.993% 
oxygen 20 ppm 
carbon dioxide 5, ppm 
hydrogen 5 ppm 
moisture 15 ppm 
argon 100 ppm 
maximum total impurities 
(non-inert) 40 ppm 


(2) Hydrogen 
Supp iweres sGanadlanmbaQuigeALl CO, ltd. 


Typical Analysis Data 


hydrogen 99n 292% 

(as remainder) 
oxygen 46 ppm 
nitrogen 31 ppm 
hydrocarbon no detected 
carbon dioxide no detected 


(Jee Nitrogen s (Low Purity.) * 
Supplier: Canadian Liquid Air Co., Ltd. 


Purity: 99.7% 


*Use for investigation of impurity effect on excess enthalpy. 


ob 
aes ui = : eo 
cae Rohe mn : 
et x08 Beal a9eno | oft % f 
| pk dee We *, 
moteyve ts port iid BA 
seer. ee | 
maq Os i 
mqag ¢ 
mag ¢ ’ gh \ 
mag cl “ 
maq Oot . pee) — 
eetsixuuqmt Letos memixem rey 
| maq O08 (tusai-non). mepestl <) 
; - : , @ 
nepoxbya” s) Vere 
‘ -b3d ,.02 1fA biupid asibened :t9biqque = 
sitsd eteylen& {sodqy? ; ‘on 
; : ; ‘oo Lo _ 
+ SRR Te meporbyd so Lm 
(tebaisme: 25) iw. ~ 
7 mqq of ; el oe i : See 
maq > ft gpotiin ee. 
Setoeseb on Pe sh eres rei 2 = ne “Jy 
bes2e75sb on sbixorb pistes <a 


*(ytitud wod) nopoxaim (8) 
-b3T-,.09 215A Biupid mekbsnsd. :teilqqae — a 


ave. naiwd ‘ an 


‘f 


47 


For various compositions the cylinders were combined as 
ratios of number of nitrogen cylinder to that of hydrogen 
Being Weto@swleto ly.ands3,. to.leform nominal 26" 5h,.wand 
76 mole percent nitrogen, respectively. 

Variations of composition with pressure for each 
combination are shown in Figure 8, which were calculated 
byseouationy (29). For this computation the compressibility 
factors were calculated from the experimental P-V-T data 
of Wiebe and Gaddy>’ for the hydrogen-nitrogen mixtures, of 
Dine for nitrogen, and of/ Michels; et al“ for hydrogen. 


Table IV shows comparisons of the composition with the 


determination by a gas chromatograph. 


Table IV 


Comparisons of Compositions Determined by Calculation and 
by Gas Chromatograph 


Puescure Mole Percent Nitrogen 
(atm) Calc. by Eq. (39) Gas Chromatograph | Devidtion 
ote 26. 0 251435 +0.65% 
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Flow Rate Calculation 


Toeuse Equation (32) for the. calculation of ‘flow 
rate, Measurements are needed for barometric pressure, P, 
room temperature, T, mixture compressibility factor at P 
and T, and time required,®, to fill out the gasometer with 
mixture, since the volume of the gasometer, Var for one 
centimeter rising of it and the rising distance, d, are 
fixed as mentioned above. 

The barometric pressure, P, and room temperature, T, 
were varied with +0.1 percent of their absolute values 
during the experiment. And mixture compressibility factor, 
Za? at standard condition is also constant, i.e., Ze 
Q29988. The error of Zn when taken as above is no larger 
than 0.15 percent when compared with Wiebe and Gaddy's 


3] Or with calculated data from the virial 


22 


experimental data 
equation of Hilsenrath, et al. 
Therefore to determine flow rate in this way, 0 is 
the only variable in the Equation (32). The flow rate 
calculated as such way, was ranged between 0.02 and 0.12 


gm-mole per second. 


ois Reproducibility of Data and Data Smoothing 


The data for the hydrogen-nitrogen system were 
reproduced with 1 percent for the entire experimental 


pressure range for each temperature level. No dependence 


E 


Geshe On the mass tlow rate was found: in this work. 
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All data were smoothed by graphical method drawing 
E , ‘ é ; : 
H’/4px(1-x) against x. The typical illustration is shown 


in Figure 9 for the temperature of =73°C. 


43° Error Analysis 


The errors involved in this data are as follows. 

The temperatures were controlled within +0.05°C for 
each temperature’ level. This can’ introduce an error 
in the excess enthalpy of 0.1 percent, which can be deter- 
mined from the Figures 16 and 17, the temperature depend- 
ence of the excess enthalpies. 

The pressure is read within 1 pound per square inch 
error which may introduce an error of 0.07 percent for high 
pressure and 0.2 percent low pressures on the excess enthalpy 
values. 

Composition determined by calculation or by analysis 
may include an error of 0.7 percent themselves which, from 
Figure 13, introduce an error on the excess enthalpies of 
1.5 percent in maximum. 

Due to the calculation of flow rate the maximum 
error on the excess enthalpies becomes 1.5 percent. 

The error on the excess enthalpies due to input 
energy calculation, was, about 0.5 percent. 

The errors due to the impurity in the pure gases, 
the time lag for the reading gauges, the pressure drop 


across the calorimeter and heat leak were negligibly small. 
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Therefore the total error introduced on the excess 


enthalpy determination from the above sources becomes less 


than 4 percent. 


5. Smoothed Data _and Comparisons with the Literature 
The smoothed data are given in Table V for each tem- 
perature level studied in each 10 atm.and 10 mole per cent 
of nitrogen. For the comparison with the literature, for 
25, 52-5, and 78 mole percent nitrogen at 20 and ee and 
forgo5 atm. at: /4 &s8,cand -73°C smoothed data were also listed. 
Comparisons with the results of Knoester, et ily are 
presented in Figures 10 and 13. Another source of the excess 
enthalpy data for this system is of the work of Hsi and The oe 
They calculated the excess enthalpies of hydrogen-nitrogen 
system from the experimental excess volume data of Zandbergen 
and Beenakker°* UusingeEquatronm e(ls/\.r higuréeésml1l - 2%and) 23 
show comparisons with Hsi and Lu's work?®, 
The agreement between this work and that of Hsi and 
Lu is within 5 percent ‘and’ with that-of Knoester,’et al. is 
Within ompenLeoentin ws (NOtCOc Then data, otekKnoester, et jal. are 


available at -72°c. Therefore they were corrected to 78°C 


for comparison with the results of this work). 


6. Data Analysis 


From a theoretical point of view, one is confronted 
with the difficult situation of presenting and describing 


ausystemeunat will cover a wide range offdensities. For 
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Figure 11. Comparison of the Excess Enthalpies 
for Hydrogen-Nitrogen System at 7°C 
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Figure 13. Comparison of the Excess Enthalpies 
for Hydrogen-Nitrogen System at -73°C 
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gases in the low density region, i.e., far removed from 

the critical state, the behaviour of the excess thermo- 

dynamic properties has been shown to obey binary colli- 

sion theory completely~. As. shown in Equation (23) the 

excess enthalpy in the region of low density is linear in 
pressure and parabolic in mole fraction at constant tem- 
perature. 

The behaviour of the heat of mixing as a function 
of pressure, for 50 mole percent nitrogen in the mixture 
with hydrogen at various temperatures studied, is shown in 
Figure 14. The excess enthalpies at temperatures of 7, 
-3,. and 73°C are shown by a linear function of pressure 
upto LPO atm .e.while at 20°C the linearity is excellent 
up to 130 atm. Figure 15 shows comparisons with the 
virial equation (Equation (23)). The excess second virial 
coefficients, E, were obtained from Brewer |. The Equation 
(23) is in good agreement at BOC, but the ‘deviation in= 
creases at low temperatures. 

Shown in Figure 13 are the isobars of the smoothed 
excess enthalpies as a function of composition at —73°Cs 
As expressed in Equation (23), all those isobars are 
excellent parabolas. 

The temperature dependence of the excess enthalpy 
is shown in Figures 16 and 17 for 50 and 80 mole percent 
WlecOgen, lLespectively, act vpressures of 20, 60, 80, and 


120 atm. with the literature values. 
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Figure 14. Isotherms of the Excess Enthalpies for 
50 Mole Percent Nitrogen of Hydrogen- 
Nitrogen Mixture as a Function of 
Pressure 
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Figure 17. Isobars of the Excess Enthalpies for 80 Mole 
Percent Nitrogen of Hydrogen-Nitrogen Mixture 
as a Function of Temperature 
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The excess enthalpy increases when the critical 
point of one of the components in the mixture is approached. 
It is interesting to observe that the excess enthalpy has 


a maximum value at the vicinity of mixture critical locus. 
B. Nitrogen-Carbon Dioxide System 


The excess enthalpies of gaseous system at the 
vicinity of the critical point of one component of the 
mixture may be determined by the liquid theories? . Also 
it is important to test the applicability of an equation of 
state and other empirical methods at that region. For this 
reason nitrogen-carbon dioxide system at 40°C was selected, 
That temperature is about 9°c above the critical temperature 
of carbon dioxide. The experiment was performed up to 130 
atm. using the same equipment as that of hydrogen-nitrogen 
system. Some additional equipment was instalied because 
of. the particular behaviour of carbon dioxide, for instance, 
Joule -Thomson effect and extreme deviation of compressibi- 
lity from ideality. The additional equipment and differ- 
ences of the experimental method from the previous system 
are described. 

Finally comparisons are made with the virial equa- 
tion of state. (Equation. (23).), the B-W-R eruetion of state, 
and Scott's liguid theories, and the Yen and Alexander's 


generalized correlation. 
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l. Material Used 

The nitrogen used for this system is the same as 
that of previous system as shown in Table III (high 
purity), and 99.9 percent pure carbon dioxide supplied 
bygcanadian= Liquid Air Co. ,,Ltd.. was used without £urther 
Dub LIacation. 

2. Experimental Work 

The commercial carbon dioxide cylinder of 50 pounds 
in net weight is in two phases (liquid and gas) at a pres- 
sure of about 68 atm. To cover the high pressure range 
it was heated up to 55°c. The experimental method was 
essentially the same as that of hydrogen-nitrogen system 
except the composition determination. 

Additional Installation 

A 55°C water bath was constructed for the high pres- 
sure carbon dioxide supply. For the random mixing of water 
in the above bath a water pump was used. 

When carbon dioxide or a mixture containing it passed 
through the regulating valve the Joule-Thomson effect was 
marked. To prevent the regulating vaive and line after the 
regulating valve from freezing due to the carbon dioxide 
behaviour, another 60°C water bath was installed between 
the operating bath (it was also called a low temperature 
bath) and the regulating Valve. Also a 6-foot-long heating 
tape wrapped the regulating valve and the copper tubing of 


the part just before the regulating valve and then the room 
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temperature bath of the previous system moved to between 
the regulating valve and the flowmeter. Since the com- 
pressibility of carbon dioxide changes strongly with 
pressure even at a constant temperature, the composition 
determination method used for the previous system was 
not suitable for this system. A gas chromatograph was 
used to determine the composition. 
Temperature Level 

The temperature studied was 40°c and it was con- 
trolled within 0.02°c during the experiment. 
Composition Preparation 

The commercial nitrogen gas cylinders were combined 
with one 50 pounds carbon dioxide cylinder to obtain various 
compositions. Figure 18 shows the combinations and the 
composition variations with pressure for various combina- 
tions. The compositions were determined by a continuous 
analyzing method with a calibrated gas chromatograph whose 
calibration data are given in Table XIV, Appendix B. 

3. Reproducibility of Data and Data Smoothing 

The experimental data of this system were collected 
as a function of pressure and at the same time the compo- 
sition variation was observed by a continuous analysing method 
as a function of pressure. Then both data were matched to ob- 
tain the smoothed excess enthalpy data. The smoothing technique 


was the same as that used for the hydrogen-nitrogen system. 
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The reproducibility of data was checked as follows. The 
data were collected for two different conditions of the 
carbon dioxide supply, one at 55°c and the other at 40°C. 
After smoothing both results, they were compared and 
were within 0.5 percent of each other. 

4, Error Analysis 

Errors introduced in the calculation of the excess 
enthalpy are from the same sources as hydrogen-nitrogen 
system, and they are almost the same magnitude. Therefore 
the error in the excess enthalpy data for this system is 
assumed less than 4 percent. 

5. Experimental Results and Comparisons 

The experimental raw data of this system are listed 
in Table XVIII, Appendix D. 

Table VI shows the smoothed data of nitrogen-carbon 
dioxide system at Oe. 

Figures 19 and 20 show isobars of the excess en- 
thalpies at 40°c and comparison with the results of 


2 Orie agreement is within 3 percent. 


Hejmadi, et al. 
At low pressures the isobars are good symmetric parabolas. 
For the increased pressures they deviate to the carbon 
dioxide side because the enthalpy change of carbon dioxide 
with pressure is a stronger function than that of nitro- 
gen. 


The pressure dependence of the excess enthalpies 


for 10, 30, 50,4,/0, and 90 mole percent of nitrogen at 
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40°C is shown in Figure 21. At sufficiently low densi- 
ties they are linear (Eq. (23)). It is obvious that the 
excess enthalpy variation with pressure is in its maxi- 
mum at around 80 atm. because the enthalpy change of 
carbon dioxide with pressure at 40°C is maximum value 
while nitrogen is changing constantly with pressure at 
same temperature. 

Figure 22 shows the excess enthalpies vanish at 
high temperatures but increase when the critical point of 


carbon dioxide is approached. 


Gee Data Interpretation and Comparisons with the 
Empirical Correlations 


The behaviour of excess enthalpy for gaseous system 
at the vicinity of the critical temperature of one compo- 
nent of the mixture, can be explained by considering the 
pressure dependence of the enthalpy for nearly critical 
fees 8 

For the nitrogen-carbon dioxide system at 40°c, 
carbon dioxide exhibits critical behaviour. The energy 
o-N, and N,-CO, 


Wi thethatLor CO,-CO. interaction. The first of these 


interactions can be neglected and we can consider the mix- 


due to the N interactions is small compared 


ing of nitrogen and carbon dioxide to its partial pressure. 


The enthalpy before mixing is then: 


x Hoo, (P,T) (40) 


-~ieneb wol 
efit t6a2 asvoiLvco 


~ixem 23i ai al 


to spaado yg 2 of3 sausosd .mis 08 Savots 3s mum 
suisyv m D +s siueesiq dtiw sbixoib nodzao 
$5 sSitteesid a instenoo pnipasdo el sepots ia siinw 
-9iutsteqmss smse 

i6 fiaiansy estaisdins easox eworle CL sispre 
70 Jtniog Jitit> soit nsdw sesetconr sud esxvssyequses Gein 
dose01qq5 ef sbhixorb modzso 

SESS ROR EISGROD DAS SOLIS IATAIOFOe Bee) oF 

enoitsloxiod Isptx tums 
Mesteye avoscEsSp 101 Briss 2e50x%S io wwotvsedsd sat 


-Oqmo>D 2 
eit pniwtsbscanos vd 


> 


[solistint 


2 
he | 0&8 ta 


Yetons 301 IuOL 


Heisgq@oo ifame ai =nr 


saon? 2160 sexrt? 
-xim sric 


<Stveesiq leistisq <i 


(08) 


; ; o 
vidnoinitive tA’ IS siypret ai nwonde ez 9°08 
at 


siveeera dittw nottsizsy yoledgns eeeoxs 


, ~ : eee 
300 10 SIUTETE 


yilasen tol 


anes | [tixD oft to yiinivtv sit 2s 
bDenmisiqnxe sd n55 ,siutxim sds te anent 
yqisditns sii to sonsbasqeb suvegetgq 
ceyeG 
; esesp 
SOrxorb 41265 -nSspoisin siz xoF 
sr {[S0iti19 a3ididxs sbixofb godwes 


aii] nofljostetat ,.00-.09 to sede ddiw 


LTSPLSHOD £62 


II ((€S) .pa&) szsentl ors yort eeis 


fI061iSI7NL .~OD—.n4 Ons 


\, 
~~ 


cM--4 srt of egb 


> 


aw bis betoslosn sd aso enotsostedak 


oJ Sbixotbh modiso bas osportin 20 pak . 


ef patxim stoied yaqlsdtas saAT 


(1.9) ght” as 
‘\ = 


a3 


30% N, 
3200 Nz - CO system 
Temperature = 40°C 
% 
2800 psi 


2400 


2000 


1600 


1200 


H® (Joules/gm-mole) 


800 


400 Vy, 


0 20 40 60 80 100 ~=—:120 
PRESSURE (atm. ) 


Figure 21. Pressure Dependence of the Excess Enthalpies 
of Nitrogen-Carbon Dioxide System 
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where x is the mole fraction of carbon dioxide. After 
mixing, the enthalpy is determined by the partial pres- 


sure of the carbon dioxide: 


x Hoos (p,T) (41) 


where p denotes the partial pressure of carbon dioxide 
in the mixture. 


Therefore, 


He (x. PAT) = x[Ho, (P/T) - Heo, (P,/T) ] (42) 
or, 

E IGS . mes 

HY (x,P,T) = x{Hog, (P+T)) - Heo, (P/T)] (43) 
where the residual enthalpy, Hace is defined as: 

H’<5 = H(P,T) - H(P=0, T) (44) 


This expression is in qualitative agreement for 
nitrogen-carbon dioxide system at 40°C, 
VirloleeguatLOMmsoOls State 

The virial equation of state truncated after second 
term is linear (Eq. (23)). A comparison is made for 50 
mole percent nitrogen with Equation (23) in Figure 23. 
To calculate the excess second virial coefficient, E, in 
Equation (23), BNo) BCO>, and Byo-CcO2 were obtained from 


Dymond and smith’, Huff and Reed*?, and Cottrell jvet Saige 


ey 


xetit&  .ebixoib nodzsp to noistost? sfom edt et * etenw 
-ae1q Isitxsq ody yd Senimesteb ei yqisdins sds , eatxia 


:sbixolb aodyas edd To eiwe 
(16) (T, a) ,008 x 


ebixolB noduso to suueesia Isit1sq sit astonsbh g sienw 
ouutxin of3 at 


,2t0Tsi8enT 


(S$) ECT.) 098 - (T,q) golle + (t, 4.x) "8 
al 
7 20 
Box aor 
(EB) { (Tv, &) aOon ~ ((2,q) eopilx = (T, 9,0) 8 
2 a” 
ol » Box, | . 
:@8 boenfteb gt , H ,yqisdins [subleet sd? etenw 
aot 


(BB). (T ,O=G)H - (T,9)H = | 


102 Jnemsexrps svidsiiisup ai ai molLeesigxs elAT 
O°OR +5 meteys sbrxotbh nodiso-aepossia 
Sis8i¢ ito sottseups Lmiaiv 

bnoose testis betsonyvit stsse to rorssiups Isizr ett 
O€ tol sbsm et nontisquon A .((€S) .pa) seanbl ak ares 
C5 stupid ar (€8) nolstsupa dtiw nepoisin tnestsq eLom 
ai ,@ ,jneioiitso> [sitiv bacose sasnxs sit stsiluolso of 
moxit benietdo sx196w sOD-cy8 . bas 1004 oot « (ES) so iteupa 
te ge ,fiesstod bas ,“Sb5on bas: tiv same odie aoa 


i 
tS 


1000 


~Ny-CO, system 
Temperature = 40°C 
50% No . 
800 Virial equation 


o This work 


600 


400 


HE (Joules/gm-mole }) 


200 


0 10 20 30 40 50 60 70 
PRESSURE (atm.) 
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Virial Equation of State 
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respectively. This equation is in agreement only up to 


SOF atm. 


B-W-R_ Equation of State 

Of the large number of equations of state, the 
B-W-R equation of state is one of the best equations for 
the high density region of gaseous systems. The compari- 
sons are made in Figure 24 and 25. The constants of the 


B-W-R equation of state were from Cullen and Kober LOr 


carbon dioxide, Bloomer and Rae and Crain and Sonntag? 
for nitrogen. They were in good agreement in the entire 
range of the experiment except in the region where the 
excess enthalpies change rapidly with pressure. 


The B-W-R equation, its constants, and the mixing 


rules for the mixture are listed in Appendix C. 


Corresponding States 
1) Ligurdviheories 


The behaviour of a binary mixture at high density 
can satisfactorily be explained by the liquid theories 
of Prigogine>° and Scott om A more general survey of the 
theories on liquid mixtures is also available from 
Rowlinson> 

Assuming that the mixture and its components obey 
the law of corresponding states, it is possible to forma 
reduced equation of state with variables reduced with 


parameters of the mixture. The excess enthalpy of a mix- 


ture was calculated by using data from reference gases 


Vv 


o¢ qu ylno insmeoxps nt ei noitsups eit? .~ylevisosqeser 
ates OE 
stet2 10 hott soupy A-W-d 

eit ,etste to enoissups to szsdmun oprsl sat 20 
tol enoizsups tesd silt to sno at stste to nottaupe A-W-G 
-Lisgmoo siiT bemortava avosesp to morpex yotensh detd edz 
e3 to etnstemoo edT .cS bois bS s1xupii at. sism eas eaos 


[ ' - ¢ 
soi el dok bas nmelluD ‘mou? sxsw ste to noissups A-W-@ 


Ef ostahoe brs aisz) bas osx bas semoold ,sbixo£b modiss 
Sxiines _— at tnesSsstpsa boop at stow yedt .nasportin sz0z 
eniz sistiw noLlpst ont mi tqsoxs jnsmizsaxs sit fo spnst 
‘Jainiesoxa: itty yibiqsx spnsdo asrqlsdzas sesoxe 

paiixim eft bas .etostenoo esi ,nottsups A-W-8 ent 
.0 xibasqgA ni bsteil sxis sauvtxim edt tO? selux 


2stex2 palbnogesz10 


soitoont Binpit (1 

yitensb dpid ts stutxim yisntd s io 1wotvsidisd efit 
seizoeds biupil st vd benisigxs sd vi trotostettse 169 
edz to ysvive [sisnsp siom A 34008 Sas O° atpoplad to 
mort sidsiisvs oais ei esrutxim biupif ao ekeniie 


|  negat tual 


yedo etnsnogme> eti bas siutxim oft tedt+ paimuees 
& miIol oF sldiacog ei ti .eetsste patbaogzertos to wel ed¢ 
figiw beoubs: esldsixsv dtiw stste to noitjsupe bsoubset 
~xim 6 to yqisritus eesoxo sdT .eizutxim sft to e193 9maxEg. sy 
| me We a 
eee6p Sonsisie1 mort sisk phiew yd botslusiso enw ext 


tc 
aan © 7 


78 


N2-CO 9 system 20% No 
3200 Temperature = 40°C 
B-W-R constants used 
N» :Bloomer and Rao® 
2800 : CO2:Cullen and Kobe ?4 
vssssnees {N2 :Crain and Sonntag? 
CO :Cullen and Kobe*4 
2400 ¢ This’ work 


De) 
oO 
(oe) 
oO 


H= (Joules/gm-mole) 
poe: 
S) Ss) 


800 


400 


PRESSURE ( atm. ) 


Figure 24. Comparison of the Excess Enthalpies of 
Nitrogen-Carbon Dioxide System with the 
B-W-R Equation of State as a Function of 
Pressure 


8y 


cM FOE mateyz ¢O)-ct 
>°Ob = siutoreamel!l 

bee zinotenes H-W-4 

‘oof bao yemoold: ooo 


P+ edoN bao nellud: sOD 


O08 
ej : 
“Hoinnoc boo noi: (My... 


2 sdor OND nelluD sOD? 


| A1OWw eZidT OOnS 
| e 
r OOOs = 
3 
| i] 
e 
L 0081 > 
hal 
| a 
. Cc 
| ° 
| oad 
1 4~OOSt = 
| w 
| Lk 
| ; 
008 
H + OOb 
| | 
eee oe 
thf 0 
OS QO! O& 06 Ob Os 8) 
iimto } JHUe2RA49 
to eaeioisisand 2zesoxd st Io GOeiteqmoD .A° sxopEa a 
oat dtiw meteve sbixo: NOLS I-Teporstiyi 
10 NOLPINGT & RE S362 Yo nolkyEuys A~-W-& 
Siveesig 
“2 
= a 


N2- CO. system 


Sian 
ae Temperature = 40°C 
3200 100atm » B-W-R constants used 
N» :Bloomer and 
2Ch Rool 
CO2 :Cullen and 
2800 Kobel 4 
N»: Crain and 
af Sonntagl 3 
CO? ; Cullen and 
Kobel 4 
2400 . —— This work 
& 
© 20 
2 00 
1 
= 
(©) 
tt 
» 1600 are 
® ae ptr 
> y : 
° i 
=> 
1200 
uw 
ae 


800 


0 20 40 60 80 100 
MOLE PERCENT NITROGEN 
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with the B-W-R Equation of State as 
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which was chosen one of the mixture components. The 
features of this theory are as follows. 

The intermolecular potential energy is assumed 
POmDeuOCetheetorm p(t je ose tiag/r) storeeach,or the pairs 
present (A-A, A-B, and B-B), where r is the distance 
between the molecular centres, 5 is a characteristic 
length, fem ecachékind@of intenacting pair and « a charac-— 
teristic energy. The components A (A-A interaction) and 
Be boBeiiceracrlonjatien cOnrorm, co the same reduced 


equation of state with the reduced variables: 


3 
DA s=Gu-—s, PSen= (45) 


where k is Boltzmann's constant. 
The the residual enthalpy for both pure components 


the same function of P* and T*, therefore, we have: 


* 
Hyver alp*krTAiis = (46) 


where N is Avogadro's number. 
The same reduced equation of state is assumed to 


be walid tor the mixture. That is: 


res 

res* Po 

H (P*¥,T*) = ——— > (477) 
un Ne 


m 
For the potential the well-known Lennard-Jones 


(6-12) potential is chosen: 
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Oo i2 Oo C= 
bap(t) = 4 ean [ - (=) (48) 


where parameters fap and o are for the interaction of the 


AB 
molecules A and B of a mixture. For unlike pairs the usual 
combination rules are assumed: 


maa lL 
Cnn = F (ay + Onn) (49) 


(50) 
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(a) Single Liquid Model 

In liquids and dense gases each molecule has at every 
moment a number of neighbours. When assuming additivity of 
forces, it seems reasonable to replace the sum of pair inter- 
action by some average interaction depending on the mole 
fractions. In the single liquid model we take an average 
<o¢(xr)> over all kinds of pairs in the system assuming a ran- 


dom mixing, which gives: 


eal Lae oX, teeta e (x) Xoo ea (x) 308 (51) 
From Equations (48) and (51), one obtains: 
<o (r)> we afx,? en (CAA) = CBB, *; 
See Sens ( (CBB) CBB)®) 
rae (CBB) - a") (52) 
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and also from Equation (48), we have: 


<p(r)> = 4 <er(SS2y7* - (S)% . (53) 


From Equations (52) and (53), we have: 


a 6 6 B.D 
— (ye ynlan + 2XXpeap%ag + Xp Bp’ pp ? 
SD Ae oc a a a (54) 
x oe fe} 12 ar Ppa eS fo) Le + x 2. fe) 12 
Ree RAC AR rae ey Na) By BE PB 
and 
Qe COND 13 2 T2¥r76 
Bee he nA AN 1-7 oR ARC ABM See BBOBE 
6r2)= af ) 5 
Xx a (e} 6 + 2x_x fe} 6 + x 2 fo} 6 
Rec AR AA A~B°AB° AB Sih ee 


(5IS))) 


Now choosing one of the two components in the mixture as 


a reference gas (denoted by subscript R), from Equation 


(45) the reduced variables are: 


‘ 
ay be By . a op3 5 ! 
<e> m Ep R 
and ¢ (56) 
T* = ne = nice TR 
€ m Ep 
and from Equations (46), (47), and (56): 
res res 
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Converting to the normal units one has: 


res Lae res 
An (x, LSet aa on He (Pp rTR) (58) 
or 
yXeSi, pom) = S82 y res (_B oe % ae T ) 
™m EE Se EN, R <e> o,3 Mo 2<ee> Mm 


(59) 


from which ce may be determined if aa 


is known at 
the indicated conditions of pressure and temperature. 


Then excess enthalpy is calculated from Equation (1). 


(b) Two Liquid Model 

Whereas in the foregoing model the molecular poten- 
tials were averaged over all the binary interactions, it 
seems obvious to apply also a model in which the inter- 
actions of A-molecules with neighbours are distinguished 
from those of B-molecules with surrounding particles. So 
in the two liguid model one considers the two average poten- 
tials <¢(r)7, and <o(x) 7, Under the same conditions as for 


A 


the single liquid model, we have: 
<o(L)>y = Xybgn(¥) + Xpoag (4) (60) 
<o(r)>p = Xp opp (4) + Xn bap (¥) x (61) 


From Equation (48), (60), and (61), we have: 
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6 6,2 

Be SOATAR AA RUAB ABI 

A 12 #2 
xe fo) ct XpeaplaR 


(62) 


12 Ze iG 


Mee (e} ar Xare re) 
A AA AA B AB AB 
rte amen | (63) 


6 3 
SOON is Ey gas 


and Ser and <o>, are similar. The thermodynamic pro- 
perties of the mixture are assumed to be the molar 
averages of the corresponding properties of these two 


liquids. This model has an expansion of enthalpy of 


mixture as: 


An = Xn <H >a + Xp <H eae (64) 
where 
oan oie) oe 
res res 
<H - (Pee, ss) a= H (—— Pay 
A rag? 4) Ep R <E>y oR? m 
€ 
— 1,) (65) 
OK 
and similarly for SH wee 


Then the excess enthalpy is obtained from: 


w pat cSte o res eeuceana res 
Hate x, {<H . Hy te + x, 1 <H B Hp a 
(c) Three Liquid Model 

Since in a binary mixture three kinds of binary 


interactions are present, it is interesting to consider 


also the mixture as a system of three fluids, each 
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characterized by one kind of pair interaction. 
The residual enthalpy of a binary system is 


then expressed as: 


(67) 


res 
where HAR 


mixed interaction, characterized by the parameter € 


is the residual enthalpy of a gas with the 


AB 


and onB as: 


(68) 


The the excess enthalpy for this model is express- 


ed as: 
Hem =e ex A - s (H + H cis (69) 


This expression shows that the mole fraction de- 
pendence of the excess enthalpy calculated with the three 
liquid model is parabolic. 

Comparisons are made in Figure 26 and Table VII. 
The two liquid model is in good agreement for this sys- 
tem when nitrogen is chosen as a reference gas. The 
enthalpy data for pure gases were taken from Hilsenrath, 
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the single liquid model the excess enthalpies calculated 
are relatively small, it is not suitable to choose carbon 


dioxide as a reference gas. 


2) Yen and Alexander's Generalized Enthalpy Departure 
Ghart 


The generalized approach for the thermodynamic 


properties prediction from the empirical data was suggest- 


ed by Lydersen, Greenkorn, and Hougen?°. 


Riexanders. modified this method. Table VII shows compari- 


Recently, Yen and 


sons of the excess enthalpies of nitrogen-carbon dioxide 
system at 40°C with other prediction method. This method 


is only in qualitative agreement. 
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Table VII 


Comparisons of the Excess Enthalpies of Nitrogen-Carbon 
Dioxide System with Calculations at 40°c for 50 Mole 
Percent Nitrogen 


Drescure Excess Enthalpies (Joules/gm-mole) 
(oer This work B-W-R® B-w-R® yea 1-Lig? AM) Totaleg, 
70 818 863 876 iikek) AMOS) 890 Ay 
80 1243 ieee 1395 TOS 444 yb 12 833 
90 23-08 2 ai, 23 On PRIME REN PA ANN 3} 1620 
100 2816 2768 2783 24360" 3158" "2947 2432 
LO 2880 2827 2843 2577) 93293 33009 2467 
120 2879 2791 2808 2660 3346 3054 2442 
Note: "4 Constant Used 
Nitrogen: Bloomer and Reon 
Carbon Dioxide: Cullen and Rober” 


Constant Used 


Nitrogen: Crain and Sonntag? 


Carbon Dioxide: Cullen and Robes? 


Single Liquid Model for Nitrogen Reference 
Two Liquid Model for Nitrogen Reference 


Single Liquid Model for Carbon Dioxide Reference 


Physical Constants. Used in Calculation 


; e/k oO ie Po Mya Zu 
Material 
a (2 


(©) 
K) (A) (°K) (atm) (l/g-mole) (-) 
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SUMMARY AND CONCLUSIONS 


A flow-type isothermal calorimeter designed for the 
study of the enthalpy changes on mixing in a gaseous 
system was constructed and tested. The experiment 
carriedout in the temperature range between -73 and 
+60°C, and for the pressures up to 130 atm. 

The hydrogen-nitrogen system was chosen to test the 
calorimeter. The excess enthalpy data of this system 
were in good agreement with the literature, within 5 
percent. 

A study has been done for the system of nitrogen-carbon 
dioxide at 40°C and UpFtOm 30 ata fon the entire com— 
position range. The data for nitrogen-carbon dioxide 
system were compared with the virial equation of state, 
the B-W-R equation of state, the generalized method, and 
the liquid theories. 

The agreement with the B-W-R equation of state is quite 
good for the entire pressure range except the region of 
rapid change of enthalpy of carbon dioxide with pressure, 
around 80 atm. The virial equation is limited to low 
densities, and the correlation of Yen and Alexander is 
in gualitative agreement. The two liquid theory is 
relatively in good agreement for nitrogen-carbon dioxide 
system, while the single and three liquid models deviate 
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To determine the mixture composition the calculation 
method may be used for the systems which are composed 
of materials whose compressibilities are in good 
linearity, and for which accurate P-V-T data are 
available. Otherwise the composition must be deter- 


mined by direct analysis. 
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RECOMMENDATIONS FOR FUTURE WORK 


An accurate automatic power supply is necessary to 
reduce experimental error and labour. 

A mass flowmeter is preferred to the present gaso- 
Meter system to decrease experimental error. 

To obtain various compositions of a system, high 
pressure regulators are recommended. 

The resistance of the heater in the calorimeter should 
be increased so that systems with large excess enthal- 
pies can be studied without using excessive currents. 
A device of gas separation and compression system is 
also recommended for the studies with expensive gases. 
The vent system for the exhaust of the mixture should 
be improved particularily when toxic gases are used. 
The tests for the virial equation of state including 
up to third term and for the Redlich-Kwong equation of 
state are recommended. 

Studies on systems of. polar and nonpolar gases at con- 


dition near the critical region should be made. 
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CALORIMETER MIXING CHAMBER DRAWING 
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APPENDIX B 


CALIBRATIONS 


00 


LOH 


Table VIII 


Calibration Data for 7555-1-B Type K-5 Potentiometer 
(Serial No. 1745006) 
6 Volts range = (02,0013) — 2muV) 


OeE6 VOLE rangen +, (Un 00m se +5 0.2 ot VV.) 


O-01L60. VOLoe range. +(0.008¢"+ O21 1V) 
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Table, IX 
Calibration Data for 8163-B Platinum Resistance Thermometer 


(Serials Now 2151205 74) 


The following values were found for the constants in 
the equation: 


Lees 


t = ———— + Slag - «21 t BG 


Ay pa 


in which t is the temperature, at the outside surface of 
the tube protecting the platinum resistor, in °c on the 
International Practical Temperature Scale of 1948 and Ry 
and Ro are the resistances of the platinum resistor at t°? 


fe) é 4 
and 0 C, respectively, measured with a continuous current 


of 2.0 milliamperes and a minimum immersion of 7 inches. 


Ro = 25.564 absolute ohms 
wsis=ted’. 010392625 
6 = 1.4924 
B = 0.1092 (t below 0°C) 
0 (t above 0°c) 


The values given are determined from measurements at 
ene Eriple pointsot water, tin point, Zine point and, the 
oxygen point. 

The error introduced by using linear interpolation 
from the data computed above equation will be less than 


Gw0001 Cc. 
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Table X 


Calibration Data for 4015-B Standard Resistor 


(Serial No. 1740278) 


The four-terminal resistance value, using the two 
top thumbscrews as potential terminals, and either the 
two side thumbscrews or the amalgamated surfaces as current 


terminals, was 


OF0999 7 sohm. 


The two-terminal value, that is the resistance be- 


tween the amalgamated surfaces of the terminals, was 


0.10002 ohm. 


Resistance Uncertainty: 0.03% 


Resistor Rating: 1 watt 


P 


re 

ot 

a i 

BIS 
Src 
. 30 
SV 
iio 


mot J 
a 
4 
t 


awSiocaednudt go? 


<dmoudst ebLe 


elsenimxss 


69s edd qneswa 


104 


Table XI 
Calibration Data for 4025-B Standard Resistor 


(Serial Nos 1737342 


The four-terminal resistance value of the above 
designated standard resistor at 25 ae. was 
9.99995 ohms. 


Resistance Uncertainty: 0.001% 


Within the interval 20° to 35°c, the change of re- 
Sistance with temperature for a manganin resistor is 


accurately expressed by the equation 


Rp = Rog [1 +a (T - 25) +8 (T - 25)7] (71) 


where Re is the resistance at T°c, R is the resistance at 


25 
250 Gy and a and g are constants. Experience shows that 
the coefficients a and Bg do not change appreciably with 


time and hence need to be determined only once. 
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Table XII 


Calibration of Moseley. 71008 Strip Chart Recorder 
(Serial No. 819-02586) 


Channel 1 Channel 2 
Actual (mV) Recording Actual (V) Recording 
0 0 0 0 
10 On 0 A 137003 
20 207.0 e 2.006 
30 30.0 3 2.990 
40 40.0 4 Sr. 9 I 
50 50%. 0 5 5. 00 
60 60.0 6 6.0L 
708 70.02 7 7201 
80 S0R02 8 S202 
90 90.02 9 9.03 
100 OOF.0 10 10.02 
PEO 109.6 IESE ee 0 2, 
120 pipe Paes} 1? T2502 
£30 12928 13 13502 
140 139:.8 14 T4202 
150 1502 ils 155.04 
160 1 6OL.L 16 16.04 
ye) 70 a0 iby) 0 4, 
180 AS UFO igs: 18.06 
190 OO AZ 19 19.08 
200 199.0 20 20, 06 
PAK 2.02.0 PAA 210.2 
220 220 Ze 22.05 
230 203 AU 23 26 0'5 
240 2407.5 24 24005 
250 25 Oxee Zo 253.05 
260 260.0 26 26.05 
270 270.0 2) 24.05 
280 2-0 On 28 28.05 
290 Z29Q75 29 29505 
300 300,..0 30 30.00 
310 SOs oul 31.05 
320 So 2s 0 a2 3205 
35.0 3 3a 0 35 350 
340 Ailes 2 34 34.05 
350 3 54co 35 355,05 
360 3G. eee), 36 365.05 
a0 3 lta S7 5/205 
380 Gish) ois: Sued a) 
390 392.0 39 39.61 


400 402.0 40 40.1 
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Table XIII 


Gasometer Calibration 


Water temperature: 18°c 
Density of water at 18°c: 0.9986 gm/cc 
Sensitivity of balance: 240% l-lb. 


The Zero point of the height of water level was taken at 
the circumference 2 cm below from the lid of hemisphere 
part of the gasometer. 


; Height of Water ; Height of Water 
mptal Weight weveliaag Total Weight “7 atreay) 
we DSaes Up Down {ibs .) Up Down 
£5 - ae ae Z 0 36.9 SW Att 
20 - -2.1 220 38.9 39 
PLE - Shera 230 40.9 Alot 
30 - 0 240 43.0 ASt 
35 ~ +1.1 250 45.0 eas 
40 - PS Si 260 Ae a7. 3 
50 - Ae 270 49.1 4953 
60 6.0 Gow 280 552 Ske 
70 8.0 8.2 290 5ae2 53.4 
80 10.0 AOL Ts 300 558 3 Bee) 
90 B2t.0 ee sug) See Sith | 
100 14.2 L4e5 320 50m 59RD 
iO 16.4 16.6 325 60.4 - 
12.0 Le. 5 L826 330 61.4 Gilze6 
130 20.4 20.6 340 63.4 63.6 
140 Z2t> 22). 6 B50 6545 65.6 
15:0 24: / 24.6 360 ORS 67.6 
160 26.8 ZGro 37.0 69.5 69.7 
EFA Ageless 28.9 380 70.6 eo 
180 30.9 Sia Ee) 390 136 PSAs 
190 52.9 3.30 395 74.6 74.6 
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Correction for the gasometer 
Eloatingecap floating cap thickness 
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AV (TaD?7h0¢ De ad 


(sas cedeayena 


= 580.Ad 


ie, (Mel OMe Xejky ON tay Gye De feley er) 


Water level changed in 9 mm when 
gasometer floating distance was 


| 65 cm. 
: =—- 9mm 5S80cc lcom 
FES ONE yaaa yada Pores Tyr 
= 8.03 cc/cm 


Correlation between volume and 


floating distance 


Choosing the starting point at 
6.1 cm and the end point at 69.7 


| sar giitiAMUTEiae | cm from the calibration tabl 
eq ET la; 2 then: 7 
"i om 5 
AV 
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seuiVe a= £2, 2k + 38.03 
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Table XIV 


Calibration: ofaGas Chromatograph 


To analyze nitrogen-carbon dioxide mixture a gas 


chromatograph was installed and calibrated as follows. 


eo conattions 


Carrier gas: Helium 
Carrier gas flow rate: 68.9 cc/min 
Column used: 60-100 mesh silica gel 


ine seinchess in thes 
imehs~babDiaghbassBrubee 


Column temperature: 107°F 


Sample gas flow rate: OF 0GNS-.C.r.H. 


2. Standard Mixture EPreparacion 


The weighing method =? was used. That is, uSing an 


evacuated cylinder of known weight, standard mixtures are 
prepared by adding and weighing each gas. The weight of 
each gas is found by difference, and the number of moles 
is computed. The mole fraction is computed from the total 


number of moles. 


3% Calibration: Data 


A total of 11 samples were used to calibrate the gas 


chromatograph. The actual data points are as follows: 
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Area Ratio Mole Fraction N 


2 
Area N2 (calculated 
SiawoEn eee from weights) 
OLO.13 QO. LoS 7 
0.2436 0.27047 
0.2795 Whe, Sieg MS 
023210 0. 35113 
0.4105 0.44431 
Opa t2 1 0.54214 
026116 0.64155 
OF 6373 On 6699S 
06725 On O62 
Oc 220 O75 290 
0.8030 0.82306 


From the data above the following equation is established. 


yae= 0214034 x" + 0585966 x (72) 


where x is mole fraction of nitrogen, and y is the area 
TaciOl. 
Equation (72) reproduces the experimental calibra- 


tion data within +0.863%. 
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N2-CO, system 


Carrier gas: Helium 
Carrier gas flow rate: 68.9 cc/min. 


Column temperature: 107°F 
Sample gas flow rate: 0.06 S.C.F.H. 


y= 0.14034 x? +0,85966x 
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APPEND Xena 


B-W-R ENTHALPY DEPARTURE FORMULA USED IN 
CALCULATION, ITS CONSTANTS USED AND 


COMBINATION RULES FOR MIXTURE 


fit £ 


(1) 


(2) 


@ tor (@ .@) (es) f= 


Enthalpy Departure Equation 


2 


fe) T=2 a7 
H - = " a 
H (BURT 2A, 4c. aa 
+ $ (2bRT - 3ayy- + = auv > 
Sy, aya =2 
+ cV T (sie er Nie | 
~yy72 
- yw 7yeTV } 
(73) 
Constants Used 
Table XV 


Constants Used for the B-W-R Equation of State 


Units: 
R= 0. 


TK) 


: fe) 
atmospheres, liters, gram-moles, K 
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SO Set iG) 


Carbon Diox 


Cullen & Ko 
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ie Re SNe Te 
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he. Bloomer & Rao® Crainies Sonntag?” 
1.053642 0.872086 
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APPENDIX E 


SAMPLE CALCULATION 


eff 
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1, Bath Temperature 


in Equation (3) 


E = 54.1862 mv 
p 
E, = 19.6426 mv 
Ms = 25.564 ohms 
and R, is obtained from Equation (71) at 22°C. as: 
R, = 9.999874 ohms 
5 Yy Benoa 602)1( 9599987 4) 


(HOG 42'6)'(25 5564) 


1.079084 


From conversion table 


Penpera buce Resistance Inverse 
Deg. ¢ Ratio Difference 

ues) On aL 2522524 

20 1.079463 Pd ayia, She) 

oe Thath 20r— (1 0079.4 63,— .1L.079084), x 252.398 


19.9043°c 


Ze Flow Rate 


Loigequat ion eG32,) 


‘ean =e F222- IC Crem 
folk = ey Sie a tel) 
P = 705 mm Hg 


OC f 
} 
STUIHI = 4 
(f[€) sorsseps at 
Vm So8f[.$c = a 
gq 
Va oS$8o.eL = 24 
emo foec.2s = A 
q 
-ap OSS ts ({T) moissupa moxt bsntssdo et a 526 
emio BY8eee.e = a” 
(AT seee. Oe) CS aBL. bc) P 
eR Braet — mee on aa care ys ah y of 
(Boc.cS) (oSha,eL) 
KBB8O0eTO.L = 
sfdss norezsvno> mort 
seisval sorsietesA SsivtsisqmeTt 
sonszeltiid OLs64a > .psd 


BSE .Ses [oz2evo.s - ef 
BeE.SeS Easevo. 1 Os 


BEE. Ses x (D80CTO.L - EabeTO.h --os a 
Soi ts dtsae 


oie wold «S 

(SE) notssupe nd 
MS\99 Sss,s = Vv 
m2 [.€9 = 5 


PH mm cdOv = gq 
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Peel 21 R0gC 

Ae te Wilt he Rays 

C) =p es omni ne Oreo) .c) Sec 
oe Eee 0 2 =a gram—mole/sec 


3a Power Supply 


Lngequarion .(3/)) 


Ey eee er OV, 
E, = Oe oe my) 
Ars W = 7.356 watts or 7.356 Joules/sec 


4. Excess Enthalpy 


In Equation (38) 


W = 7.356 Joules/sec 
Di = 0.0795 gram-mole/sec 
ee egos Joules/gm-mole 


007915 


Note: 
This result is listed in the eighth line for the 
Eemperacure Of 19.93 + 0.05-C, nominal 26 mole percent 


nitrogen in Table XVII, Appendix D. 
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